A Mechanism for Hadron Molecule Production in pp(p) Collisions 
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We propose a mechanism allowing the formation of loosely bound molecules of charmed mesons 
in high energy proton-(anti)proton collisions. 
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Introduction. The problem of understanding the loosely bound hadron molecule formation in pp(p) collisions 
at Tevatron and LHC energies is still open. A recent measurement by the CMS Collaboration |1] basically 
confirms, at higher energies, older Tevatron results on the prompt production of A(3872) which were first 
addressed in Q. Looking at these new results the questions remain the same as those raised in @|: how is 
that possible that a very long lived molecule of a D° and a D*° meson, with binding energy compatible with 
zero, could be formed within the bulk of the hadrons ejected in very high energy pp(p) collisions? Is it the 
A (3872) that molecule? 

The reply given in Q to the former question was sharply negative. In that paper we performed numerical 
simulations with standard hadronization algorithms (Herwig and Pythia) tuned to fit data on the production of 
open charm mesons and sought D°D*° pairs with reasonably low relative momentum in their centre of mass so 
as to be eligible candidates for becoming molecular loosely bound states. The number of selected pairs allowed 
to estimate an upper bound on the prompt [l2| production cross section of the A(3872) which was found to be 
at least 30 times smaller than the experimental value. 

Our analysis was reproduced, with similar results, in Q, where it was also observed that a more appropriate 
treatment of Tevatron data would rather indicate a discrepancy with theoretical expectations by a factor of 300. 

Such a gap did not seem to be unbridgeable to the authors of Q , who resorted to final state interaction (FSI) 
mechanisms in the D°D*° system in order to improve the theoretical cross section up to the experimental value. 
The approach there used was criticised in Q leaving the controversy somewhat unsolved 

Molecular A(3872). On the other hand, during the last few years, the idea of a molecular X, in diverse 
incarnations @, has been corroborated by the lack of observation of its nearly degenerate charged partners, 
required by the antagonist tetraquark model [tJ ■ For these reasons we come back here to the problem of the X 
formation in high energy hadron collisions being motivated by a completely different approach. 




FIG. 1: The elastic scattering of a D° (or D*°) with a pion among those produced in hadronization could 
reduce the relative momentum fco in the centre of mass of the D°D°* pair. 



In our view the X could rather be the meson-molecule analogue of the stable deuterium. 

Given the large number of pions produced in the neighbourhood of the open charm meson pairs in momentum 
phase space, it is plausible that some of those pions could scatter elastically on the D° or D*° component of 
the would-be- molecule changing the relative momentum in the centre of mass of the pair, kg, towards lower 
values - see Fig. Q] We can assume the initial total energy £ of the pair to be positive. However, if fco = |fco| 
gets smaller due to an interaction with the pion, £ might be found shifted down to some negative - close to 
zero - value, provided that the D°D*° pair is under the influence of some (unknown) attractive potential, say 
a square well potential, similar to the simplest description of deuterium. 

In these respects the X would be a genuine, negative energy, bound state of D°D*° whose lifetime is en- 
tirely regulated by the lifetime of the shorter lived component D*°; we would estimate then a total width 
rtot(A) ~ 65 keV 0. There are no energetic arguments to stabilize the D* in the attractive potential. 



Such a mechanism is therefore somewhat opposite to the one based on FSI, where the D°D*° pair should 
rescatter remaining isolated from other hadrons potentially produced close in phase space [H, 0| • 

One more reason to pursue the approach described above is that the resonant scattering D°D*° — >• X — >• 
D°D*° is difficult to be reconciled with the general expectations that can be drawn for the total scattering cross 
section of two particles allowing a shallow bound state with energy \£\ ~ 0, as described in the 'Low equation' 
formalism, see fsj]. Resonant scattering D°D*° — >• X — >• D°D*° can be computed using available data on X 
decay branching fractions (in order to compute the g XDD , coupling) and averaging the cross section, (a) /, with 
the distribution /(fco) of D°D*° pairs obtained by hadronization algorithms [13| |. It is only when ko is smaller 
than some critical value A that the resonant scattering into X has a non negligible probability to occur. We 
find a scattering legth a of about 4 fm for a total width T(X) ~ 1.2 MeV (the g XDD , coupling is a function of 
the X total width) and A = 50 MeV. The scattering length a decreases for smaller values of the total width - 
see Fig. [5] Shifting A towards higher values, A ~ 10 GeV, a decreases to few cents of a fermi. 

On t he oth er hand, the scattering length expected for scattering with a shallow bound state is 
a = H/ y/2[i\£ | ~ 12 fm {£ ~ —0.14 MeV). Such a result, as discussed in is independent on the (unknown) 
scattering potential. 
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FIG. 2: Scattering length for the D°D*° -> X -> D°D*° process as a function of the X total width. The initial pairs 
are selected with |feo| < 50 MeV which, in our simulations, represent a few parts over 10 5 of the total. The error bands 

account for uncertainties on the data we used. 



Analysis method. The binding energy of the X is estimated from the mass difference with its constituents 
£ ~ —0.14 ±0.22 MeV. A discrete level at this energy (take the central value) can be accommodated in a square 
well with a depth of about —7 MeV [lH and a range r$ ~ 3 fm. 

Let ipo be the wave function ass ociated to this level. The average size of the molecule is found to be 
v/ (r 2 )^ ~ 10 fm and a value of y/ (feo)v'o — 50 MeV is determined. Those pions scattering elastically on 
D° or D*° and making the fco of the pair lower than 50 MeV are able to drop the total energy down to £ 
and form a genuine D°D*° bound state. It is our purpose here to seek such pions and to study numerically 
their elastic interactions with the D° or D*° mesons adapting standard hadronization tools such as Herwig and 
Pythia. 

As discussed first in Q, the spectrum of D°D*° pairs can be represented by a monotonically rising histogram 
in fco. Because of the interaction with pions, pairs with high relative COM (centre of mass) momenta, the 
majority, could either be pushed to higher momenta or to lower ones. If even a small part of them were 
rearranged within lower relative momenta, there could be a significant effect of feed-down of pairs towards lower 
bins, even in the far low energy region below 50 MeV. Populating that region means increasing the formation 
probability of the loosely bound X. 

To perform a first qualitative exploration of this phenomenon, we start by generating samples of pp — > cc 
events in Herwig and Pythia, at Tevatron COM energies (y/s = 1.96 TeV). We list the events containing D°D°* 
(resp. D°D *) as a function of fco- The cuts imposed at parton level are: pj_ > 2 GeV and |y part | < 6. 

The distributions da/d(Aip), where Aip is the difference in azimuthal angles between D° and D*~ , as discussed 
in 0, are reproduced by choosing the following cuts on the final mesons: open charm meson pairs have 
5.5 GeV < p± < 20 GeV and \y\ < 1. These cuts allow to reproduce very well CDF data on da/d(Aip) if 
a full Quantum Chromodynamics (QCD) generation of events is performed (cc + gg + gq + qq...). D°D*° pairs 
in the bin A<p = [0°, 18°] are the main would-be-molecule candidates. We observe here that the numerical 
generation of pp —> cc partially fills (« 10%) the = [0°,18°] bin with respect to the full QCD one. In 
addition, in the central region, which is enforced by the cuts, we have to match our results with those of some 
Matrix Element Monte Carlo, like Alpgen [l(J, more than just using shower algorithms. We will present the 
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results of the full QCD simulation, which is much more time consuming, in a future paper. 

To optimize the selection of events, we choose the 10 most complanar pions to the D°D*° plane, then we 
randomly choose the meson the pion will interact with (say the D°), and finally we select the most parallel pion 
to the non-interacting meson (say the D*°) - see Fig. [TJ In physical events, we expect such a pion to be the 
most effective one to the phenomenon we are describing. 

The elastic interactions with the pions are regulated in the ttD° COM by the matrix elements 

(n(p)D(q)\D*(P,r 1 ))=g^ DtV -p 
(K(p)D*(qA)\D*(P,r))) = ^gf e a ^ s X a r,^q 5 

where the couplings used are g itDD , ~ 11, g mD , D * ~ 17, see [ll|- After the interaction with the pion has taken 
place in the COM D°tt frame, we boost back the D° in the laboratory (LAB) frame and check if the 'new' 
D°D*° pair passes the cuts we fixed for the final meson pairs. 

We can trace, event by event, the variation fco k' Q of each D°D*° pair filling a 2D histogram of transition 
probabilities V (ko, Afco). Since the interaction with pions can change the p± and y of the molecule, a pair 
might fail the strict meson cuts before the interaction and pass them after it (a 'gained' would-be-molecule) 
and viceversa (a 'lost' one): see Fig. [3J 

The open charm mesons might interact with pions more than once before a molecule is formed. Roughly 
speaking the ttD° — > ttD° scattering is proportional to g^-riD* whereas the D* — > Dir decay is 'slower' by 
9kDD* G3- We assume that a single D°(D*°) might scatter, on average, with 2-^3 pions before the relative 
distances among the flying-out hadrons are such that the interactions are suppressed [16j. 

Therefore, for each pair, we wish to evaluate k^ after n interactions. We do it according to the probability 
distribution functions (PDF) as extracted from V (fco, Afco). We build a set of PDFs Vi (Afco) for each bin i in 
da/dko. We assume that the PDFs will be the same for all the interactions, like in a Markov chain. For each 
event we have a k$ , falling in some particular bin i. We randomly extract a Afco according to the distribution 

Vi (Afco) and sum |fc ™' ) + Afco| = fcg™ +1 ^ thus producing a new histogram. 

We must also take into account the 'lost' and 'gained' would-be-molecules. In each iteration, we generate the 
number of 'lost' and 'gained' ones, l^ n \ g^ n \ according to Poissonian distributions with mean values g^\ 
We implement the following algorithm: i) before the n-th interaction, we drop out a number of pairs, ii) 
we produce the new histogram as a result of the interaction with one more pion, Hi) after that, we decide to 
'gain' a g( n > number of pairs. 

Results. The results are showed in Fig. [3J The bin we are more interested in is the first one, with 
fco < 50 MeV. The number of pairs obtained for that bin are reported in Table I. As one can see from these 




FIG. 3: Number of D°D*° pairs (events) counted with Herwig (left panel) and Pythia (right panel) when generating 
10 10 pp — > cc events at y/s = 1.96 TeV with the cuts on partons and hadrons described in the text. The Ow histogram 
reproduces the shape found in Q| . The histograms named Itt and Sir are related to the elastic scattering of open charm 
mesons with one or three pions selected as described above. In the insects we report a broader fco range. 



plots the feed-down mechanism towards lower relative momentum bins is very effective once the interaction of 
a D or a D*° with a pion from the hadronization is taken into account. The effect gets magnified if successive 
interactions are allowed (up to three). In the insects we show a broader range in fco. It is evident here that the 
elastic scattering with a pion is also causing a net increase of would-be-molecule pairs: it forces a number of 
pairs to pass the p± > 5 GeV and \y\ < 0.6 cuts, which otherwise would be failed. 
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k < 50 MeV 


Ott 


In 


3tt 


Herwig 


10 


19 


802 


Pythia 


3 


21 


814 



TABLE I: The population of the the k < 50 MeV bin (D°D*° pairs), after 0, 1, 3, tt£>M interactions. 

The results showed in Table U are indicating qualitatively that the mechanism described in this letter indeed 
occurs in numerical simulations of pp collisions and might play an important role in physical events. For a full 
determination of prompt production cross sections we need to switch from pp — > cc to the full QCD generation 
pp — > cc + gg + gq + qq... which is a harder task in terms of numerical computation, yet, from the exploration 
here reported, we have a clear clue on what to expect. 

Conclusions. We have presented a new mechanism to explain the prompt formation of loosely bound open 
charm meson molecules at hadron colliders as induced by elastic scattering with comoving pions. Simplified 
numerical simulations show that pions produced in hadronization might be effective at decresing the relative 
momentum in the center of mass of the D meson pair which, if under the influence of an attractive potential, 
might therefore be found at some small negative energy, like in a shallow bound state in a potential well. Such a 
bound state will have a lifetime which is as long as the D*° one, T ~ 65 keV, still well below actual experimental 
resolution. With the results of the full numerical simulations we will provide expected prompt cross sections 
for the production of the X(3872) at the LHC. 

Considering the known limits of the available hadronization models, the results of numerical simulations have 
to be taken as compelling but qualitative descriptions of the suggested mechanism. We believe that several 
more investigations in this direction are possible. 
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